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ABSTRACT 


Direct  heat  addition  of  flowing  gases  by  radiation  absorption  is 
considered.  In  the  visible  wavelength  regime,  the  interaction  between 
concentrated  solar  energy  and  a  flowing  gas  is  being  modeled  numerically. 
Additional  analyses  of  1  aser-gasdynami c  interactions  are  also  being 
considered.  Implicit  time-iterative  procedures  originally  developed  for 
transonic  flows  are  being  adapted  to  the  low  Mach  number,  low  Reynolds 
number  regimes  of  interest.  Two-dimensional  solutions  show  that  absorpti 
plasmas  can  exist  over  a  wide  range  of  flow  speeds  and  that  buoyancy  is  a 
dominant  factor  in  the  forced  convection  flowfields.  In  the  microwave 
regime,  an  experimental  investigation  of  various  absorption  modes  is 
underway.  An  overview  of  the  microwave  resonant  cavity  experiment  and 
initial  experimental  results  with  nitrogen  is  presented.  Numerical 
predictions  of  the  planar  and  propagating  plasma  mode  of  microwave 
absorption  in  hydrogen,  helium  and  nitrogen  are  also  discussed. 


Statement  of  Work 


Determine  experimentally  the  conditions  under  which  plasmas  can  be 
initiated  and  sustained  in  the  three  energy  addition  modes  (TM^,  TEyj, 
and  planar).  Initial  testing  will  be  with  nitrogen  and  helium. 

Nitrogen  will  be  used  to  simulate  the  molecular  nature  of  hydrogen 
while  helium  will  simulate  atomic  hydrogen.  Final  testing  will  be  with 
hydrogen.  Parameters  to  be  examined  include  gas  composition,  pressures 
and  flow  rates  and  microwave  power.  Quantify  system  heat  losses. 
Measure  spectroscopically  electron  and  ion  temperatures  and  densities 
in  microwave  generated  plasmas  in  the  three  energy  addition  modes.  Due 
to  high  gas  pressures,  local  thermodynamic  equilibrium  (LTE)  will  be 
initially  assumed  although  nonequilibrium  effects  will  be  examined. 

The  principal  measurement  will  be  ion  temperature  as  this  translates  to 
thrust.  The  effects  of  various  hot  and  cold  gas  mixing  schemes  on 
final  temperature  will  be  studied. 


Status  of  Research  F.ffort 


Free-floating  spherical  plasmas  have  been  generated  in  stationary 
nitrogen  gas  contained  inside  a  spherical  quartz  vessel  located  within  a 
cylindrical  resonant  cavity  operated  in  the  TMq^  mode.  The  plasmas  are 
approximately  two  inches  in  diameter  and  are  centered  within  the  five  inch 
ID  quartz  sphere.  Power  coupled  to  the  plasma  was  measured  as  a  function 
of  nitrogen  pressure  and  microwave  power  input  to  the  cavity  only  for  those 
combinations  of  pressure  and  power  which  resulted  in  the  plasma  being 
stabilized  in  the  center  of  the  quartz  sphere  away  from  adjacent  walls. 
Figure  1  plots  the  coupling  efficiency  as  a  function  of  absolute  pressure 
for  several  input  power  levels.  The  coupling  efficiency  is  defined  to  be 
the  percentage  of  input  microwave  power  actually  absorbed  in  the  gas.  It 
can  be  seen  that  the  coupling  efficiency  rises  with  increased  pressure  but 
decreases  with  increased  power.  Figure  2  plots  the  total  microwave  power 
absorbed  by  the  gas  and  shows  that  the  total  power  absorbed  initially 
increases  with  both  increased  input  power  and  gas  pressure  but  levels  off 
above  input  powers  of  350  W  and  gas  pressures  of  15  kPa.  Figure  3  plots 
the  maximum  coupling  efficiency  obtained  from  Figure  1  as  a  function  of 
input  power  and  shows  the  coupling  efficiency  dropping  from  80%  to  30%  as 
the  input  power  is  raised  from  100  to  500  W. 

A  numerical  model  of  the  one-dimensional  planar  propagating  microwave 
plasma  in  hydrogen,  helium  and  nitrogen  gas  has  been  successfully 
formulated. 1  This  model  numerically  integrates  the  system  of  governing 
equations  consisting  of  the  one-dimensiona l_ steady  energy  equation 


pu  C 
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dT 

dx 


(1) 
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where  p  is  Che  gas  density,  u  is  Che  gas  velocity,  Cp  is  the  specific  heat, 
T  is  the  gas  temperature,  A.  is  the  thermal  conductivity,  0  is  the 
electrical  conductivity  and  E  is  the  electric  field  vector;  and  Maxwell's 
equation  describing  the  propagation  of  the  microwave  energy 


dE 
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where  e  is  the  permitivity  and  p  is  the  permeability.  The  specific  heat, 
thermal  conductivity  and  electrical  conductivity  are  functions  of 
temperature.  The  electromagnetic  energy  is  absorbed  by  the  plasma  as  the 
temperature  rises  with  some  microwave  energy  being  reflected  or  transmitted. 
Due  to  thermal  conduction  to  the  cold  gas  ahead  of  it,  the  plasma 
propagates  toward  the  microwave  energy  source  at  a  velocity  determined  by 
the  energy  balance  between  the  absorbed  microwave  power  and  the  heated  gas 
which  is  convected  away  downstream. 

The  two  governing  equations  were  numerically  integrated  using  a 
fifth/sixth  order  variable  step  Runge-Kutta  scheme.  An  iterative  method 
was  used  to  determine  the  propagation  velocity  eigenvalue,  pu,  similar  to 
the  method  used  by  Kemp  and  Root^  to  solve  for  the  propagation  velocity  of 
the  laser  heated  plasma.  The  propagation  velocity,  maximum  temperature  and 
percent  power  absorbed  were  calculated  as  functions  of. the  input  microwave 
power.  Figure  U  plots  the  propagation  velocity  as  a  function  of  input 
power  and  shows  the  propagation  velocity  rising  with  increased  microwave 
power.  The  decreased  velocity  for  the  higher  gas  pressure  is  primarily  due 
to  the  increased  gas  density  since  pu  is  the  propagation  eigenvalue  and  is 
almost  equal  for  the  two  gas  pressures.  Figure  5  plots  the  maximum  gas 
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temperature  as  a  function  of  input  power  while  Figure  6  plots  both  the 
percent  power  absorbed  and  the  percent  power  reflected  by  the  propagating 
plasma  as  a  function  of  input  power.  The  maximum  gas  temperature  is 
constant  for  hydrogen  and  nitrogen  because  significant  dissociation  is 
occurring  in  this  temperature  range  and  is  absorbing  the  additional 
absorbed  power.  The  maximum  gas  temperature  for  helium,  which  is  not 
dissociating,  is  seen  to  rise  with  increasing  input  power.  It  was  found 
that  all  of  the  input  power  was  either  reflected  or  absorbed  with  no  power 
being  transmitted  through  the  plasma  and  that  the  percent  power  absorbed 
for  hydrogen  decreases  with  increased  input  power  while  the  percent  power 
absorbed  for  helium  and  nitrogen  remained  fairly  constant  as  a  function  of 
input  microwave  power.  The  results  of  the  numerical  model  will  be  compared 
to  experimentally  measured  values  as  they  become  available. 
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Coupling  Efficiency  Versus  Absolute  Pressure  for  Various  input  Powers 


Fig.  X  Coupling  efficiency  of  a  free-floating  spherical  plasma  in  stationary 
nitrogen  gas  as  a  function  of  absolute  pressure  for  several  input 


owe r  Coupled  to  the  Plasma  Versu 


Fig.  2  Total  microwave  power  coupled  to  a  free-floating  spherical  plasma  in 
stationary  nitrogen  gas  as  a  function  of  absolute  pressure  for 
several  input  power  levels. 
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I.  STATUS  OF  RESEARCH  EFFORT 


A.  Computational  Algorithm 

The  advanced  computational  algorithms  that  have  been  developed  for 
transonic,  external  aerodynamics  in  the  past  decade  are  being  adapted  for 
low  speed  flows  with  strong  heat  addition.  The  specific  area  of  focus  has 
been  on  the  computation  of  the  interaction  between  concentrated  solar 
radiation  and  flowing  gases,  but  the  techniques  are  equally  applicable  to 
laser-gas dynamic  interactions  and  microwave  absorption,  as  well  as  to  the 
more  traditional  combustion  problems  that  occur  in  propulsion  environments. 
The  family  of  techniques  being  addressed  is  time  dependent  procedures. 

These  are  attractive  because  they  apply  to  either  viscous  or  inviscid  flows 
and  provide  high  accuracy  and  efficiency.  In  addition,  time  dependent 
procedures  are  applicable  to  either  steady  or  unsteady  flows,  making  them 
candidates  for  studying  both  the  steady  flow  characteristics  of  propulsion 
environments  as  well  as  the  linear  and  nonlinear  stability  characteristics 
of  these  flowfields. 

The  extension  of  time  dependent  algorithms  from  external  transonic 
flows  to  internal  propulsion  environments  required  advances  in  several 
directions.  The  methods  had  to  be  modified  to  handle  very  low  Mach  number 
flows  (which  are  nonetheless  compressible  because  of  the  energy  deposition) 
in  addition  to  the  transonic  flows  for  which  they  were  developed.  In 
addition,  they  had  to  be  extended  from  high  Reynolds  number  applications  to 
low  Reynolds  number  conditions  because  some  of  the  energy  deposition 
processes  described  above  take  place  at  low  Reynolds  numbers.  The 
modifications  which  enabled  these  low  Mach  number  and  low  Reynolds  number 
extensions  are  discussed  below. 
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The  equations  of  motion  for  a  flowing  gas  with  heat  addition  (by 
absorption  of  radiation,  I^R  heating,  or  combustion)  can  be  expressed  in 
vector  form  as. 
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Here,  Q,  E,  F,  H,  Ev  and  Fv  are  vectors.  The  primary  dependent  variables 
are  contained  in  Q  while  the  remaining  five  vectors  are  nonlinear 
functions  of  Q.  These  vectors  are  defined  as, 
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where  Qv  *  (P»  u,  v,  T)T,  and  Rxx,  RXy,  RyX,  and  Ryy  are  matrices 
containing  only  the  viscosity  and  the  thermal  conductivity.  For 
computational  purposes,  the  equations  are  transformed  to  generalized 
nonorthogonal  coordinates.  Convergence  to  a  steady  state  corresponds  to 
following  Eqn.  1  through  a  transient  from  some  arbitrary  initial  condition 
to  a  final  steady  condition.  We  begin  by  discussing  the  inviscid  terms  on 
the  left-hand  side  and  then  add  the  viscous  terms  on  the  right-hand  side. 

An  example  of  the  manner  in  which  typical  time  dependent  methods  are 
slowed  down  as  the  flow  Mach  number  is  reduced  is  given  in  Fig.  1.  The 
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reason  for  the  decreased  efficiency  is  because  the  eigenvalues  of  the 
system  become  increasingly  stiff  as  Mach  number  is  reduced.  We  have 
developed  two  methods  for  circumventing  this  problem.  In  the  first,  the 
time  derivatives  are  multiplied  by  some  matrix  T  which  is  chosen  in  such  a 
way  that  the  eigenvalues  remain  well  conditioned  at  low  Mach  numbers.  This 
matrix  preconditioning  procedure  corresponds  to  modifying  the  physical  time 
derivatives  to  control  the  acoustic  speeds  and  provide  faster  convergence. 
The  convergence  of  an  inviscid  problem  with  this  preconditioning  is  shown 
in  Fig.  2  for  Mach  numbers  down  to  10“5.  As  can  be  seen  from  the  figure, 
preconditioning  allows  us  to  maintain  rapid  convergence  to  low  Mach  numbers 
but  eventually  round-off  errors  are  encountered  in  the  pressure 
computations.  These  round-off  errors  prevent  the  use  of  the  procedure  at 
Mach  numbers  below  about  10~4.  A  summary  of  some  initial  preconditioning 
results  is  given  in  Ref.  9. 

The  second  method  for  accomplishing  low  Mach  number  computations  has 
been  to  use  a  perturbation  expansion  of  the  equations  of  motion  to  obtain  a 
low  Mach  number  system.  In  this  low  Mach  number  system,  artificial 
acoustic  modes  are  introduced  to  provide  a  well-conditioned  eigenvalue 
system  at  any  Mach  number.  This  system  also  includes  a  rescaling  of  the 
pressure  and  all  inviscid  applications  have  shown  it  to  be  very  effective. 
We  have  again  obtained  Mach  number  independent  convergence  to  Mach  numbers 
as  low  as  10"6  which  exceeds  the  low  speed  range  of  the  laser  propulsion 
problems  of  interest  and  the  pressure  scaling  eliminates  the  premature 
round-off  errors  noted  with  preconditioning.  The  low  Mach  number  expansion 
procedure  is  described  in  Refs.  9,  10  and  12.  Representative  convergence 
rates  for  the  low  Mach  number  scheme  are  given  in  Fig.  3. 
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The  addition  of  viscous  terms  to  this  low  Mach  number  algorithm 
results  in  slightly  enhanced  convergence  rates  at  higher  Reynolds  number 
(Re  -  1000)  as  compared  with  inviscid  conditions  (Re  -  ®)  but  at  lower 
Reynolds  numbers  the  convergence  rates  begin  to  decrease  again  as  shown  on 
Fig.  4.  Unfortunately ,  Re  -  10  is  representati ve  of  laser  absorption 
conditions.  These  reduced  convergence  rates  at  low  Reynolds  numbers  can  be 
avoided  by  proper  rescaling  of  the  time  step  (as  discussed  in  detail  in 
Ref.  12).  Representative  results  for  an  optimum  time  step  are  given  in 
Fig.  5.  The  convergence  rates  given  here  are  for  Mach  number  and  Reynolds 
number  conditions  representati ve  of  those  in  solar  or  laser  absorption 
regimes . 

An  example  of  flowfields  calculated  with  this  algorithm  for  Reynolds 
numbers  from  infinity  (inviscid)  to  Re  -  0.05  are  shown  on  Fig.  6  for  flow 
in  a  duct  with  a  specified  rate  of  volumetric  heat  addition.  These  various 
flow  regimes  include  all  conditions  expected  in  beamed  energy  propulsion 
systems  from  the  energy  absorption  regime  through  the  throat.  Similar 
calculations  for  energy  addition  through  wall  heating  are  given  in  Fig.  7. 
The  overall  domain  of  convergence  in  a  Mach  number,  Reynolds  number  plane 
is  shown  on  Fig.  8. 

B.  Application  to  Laser  Absorption  Problem 

Applications  of  the  modified  algorithm  to  the  laser  absorption  problem 
are  given  in  Figs.  9-13.  Figure  9  shows  a  representative  grid  with 
moderate  stretching  and  a  superimposed  ray-traced  laser  beam. 

Representative  convergence  rates  for  the  laser  absorption  calculation  are 
shown  in  Fig.  10  for  three  cases:  a  direct  method  in  which  the  laser  rays 
and  the  fluid  dynamics  equations  are  solved  directly  and  simultaneously; 
one  with  direct  inversion  of  the  fluids  equations  in  which  the  laser  ray 
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calculations  are  uncoupled;  and  one  which  uses  the  AOI  method  for  the 

i 

f 

fluids  equations  and  again  is  uncoupled  from  the  laser  beam.  Clearly,  the 
laser  ray  calculations  should  be  coupled  with  the  fluids  equations,  but  the 
CPU  requirements  for  the  complete  system  remain  quite  large  even  though 
convergence  is  achieved  in  three  to  five  iterations.  Details  of  these 
techniques  are  to  appear  in  Ref.  13. 

Figure  11  shows  temperature  contours  for  laser  absorption  in  argon  at 
several  laser  power  levels.  These  calculations  are  for  a  flow  velocity  of 
1  cm/s,  and  ignore  the  effects  of  gravity.  The  effects  of  varying  the 
inlet  velocity  at  fixed  laser  power  are  seen  to  result  in  much  greater 
movement  of  the  heating  zone  as  demonstrated  in  Fig.  12,  while  Fig.  13 
demonstrates  the  effect  of  scaling  up  the  laser  size.  This  latter 
variation  does  not  show  strong  movement  of  the  plasma  region,  but  size 
-  scaling  beyond  that  indicated  appears  to  result  in  conditions  which  no 
longer  allow  a  steady  solution.  Present  emphasis  here  is  on  improving  the 
radiation  loss  model  and  Investigating  the  effects  of  size  scaling  in  more 
detail.  The  radiation  loss  model  is  described  below. 

C.  Absorptivity  Characteristics  of  Hydrogen  Alkali  Mixtures 

To  enable  absorption  predictions  of  solar  energy,  we  are  developing  a 
broadband  model  of  the  absorption  characteristics  of  hydrogen-alkali 
metal  vapor  mixtures.  Alkali  metal  vapors  are  known  to  be  efficient 
broadband  absorbers  of  visible  wavelengths.  Theoretical  studies  have  shown 
that  the  absorption  of  concentrated  solar  energy  by  alkali  metal  vapors  is 
due  to  two  mechanisms:  photoionization  and  dimer  transitions.  The  dimer 
absorption  bands  of  metal  vapors  extend  over  most  of  the  solar  spectrum. 

For  sodium  and  potassium,  the  strongest  bands  lie  near  the  peak  of  the 
solar  spectrum  while  for  cesium  these  bands  extend  into  the  infra-red 
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region.  When  mixed  with  hydrogen,  the  alkali  metal  vapors  will  form 
hydrides  that  are  also  effective  absorbers  in  the  visible  regime,  but  these 
additional  absorption  mechanisms  are  being  ignored  for  now. 

As  an  approximate  representation  of  the  absorption  characteristics  of 
alkali  metals,  we  are  attempting  to  model  the  radiative  transfer  in  various 
wavenumber  regimes  as  a  series  of  bands.  In  each  band,  the  molecular 
absorption  coefficient  is  based  upon  the  Franck  Condon  principle  and  uses  a 
method  called  the  quasi-static  line-broadening  principle.  (From  a 
classical  viewpoint,  the  Franck  Condon  principle  states  that  spontaneous 
transitions  between  electronic  states  are  most  likely  to  occur  when  the  two 
nuclei  are  at  their  extreme  vibrational  positions.  Further,  it  assumes 
that  the  accompanying  photon  absorption/emission  takes  place  so  rapidly 
that  the  positions  of  the  nuclei  do  not  change.)  This  quasi-static 
approximation  requires  the  knowledge  of  the  potential  energy  curves  for  the 
different  molecular  states  and  the  oscillator  strengths  corresponding  to 
each  transition. 

Under  the  quasi -static  approximation,  the  probability  of  a  transition 
occurring  at  frequency,  v,  while  the  nuclei  are  separated  by  a  distance,  R, 
(here,  hv  is  the  energy  difference  between  the  upper  and  lower  states), 
can  be  used  to  estimate  the  molecular  absorption  coefficient: 


where  f  is  the  oscillator  strength,  g  is  the  multiplicity  of  the  lower 
electronic  state,  V(R)  is  the  binding  energy,  and  other  symbols  have  their 
standard  meaning.  The  derivative  dR/dv  in  this  expression  is  computed  from 


the  potential  energy  diagrams  of  the  molecule.  Calculations  based  on  Eqn. 
2  for  sodium  shown  good  agreement  with  experimental  spectra. 

Although  we  are  planning  to  use  this  approach  as  an  interim,  we  are 
using  an  alternative  procedure  based  upon  a  fit  to  existing  experimental 
data, 

2  2 

yv.T)  -  Hr  *9  N2  G(m)  exp  C-^i]  (3) 

where  G(v)  is  the  line  shape  factor.  For  quasi-static  broadening,  this 
profile  is  torentzian. 
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Here,  Av  -  v  -  v0  and  y  is  the  linewidth  parameter  which  is  taken  from 
experimental  data.  Representative  absorption  coefficients  for  sodium  and 
potassium  are  given  in  Figs.  14  and  15  for  various  temperatures. 

For  the  range  of  temperatures  of  interest  (3000K  to  5000K)  broadband 
absorption  of  solar  radiation  by  bound-free  or  photoionization  transitions 
may  also  be  important  in  the  visible  region.  The  semi-classical  derivation 
for  the  bound-free  absorption  cross-section  yields  the  Kramer’s  formula. 
Combining  this  with  the  Boltzman  equation  for  the  number  of  atoms  in  the 
nth  state,  we  obtain  the  bound-free  absorption  coefficient  as, 
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where  E]  is  the  energy  of  the  electron  in  the  ground  state  and  En  is  the 
energy  of  the  n^*1  state.  Figure  16  shows  the  bound-free  absorption 
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coefficient  of  Na  at  a  temperature  of  5000K.  The  total  absorption 
coefficient  of  an  alkali  metal  vapor  is  the  sum  of  the  dimer  and  bound-free 
coeffi cients . 

0.  Radiation  Loss  Modeling 

Once  the  absorptivity  characteristi cs  of  the  working  fluid  are 
determined,  the  next  step  is  to  develop  a  model  for  predicting  the 
radiative  losses  from  the  hot  zone.  Because  the  radiative  flux  that 
appears  in  the  energy  equation  arises  because  of  a  volumetric,  rather  than 
a  local,  effect  its  evaluation  requires  inordinately  large  amounts  of 
computer  time  and,  consequently,  it  must  be  modeled  rather  than  computed. 
This  volumetric  dependence  also  means  that  no  simple  scaling  laws  can  be 
developed  for  the  general  case. 

The  engineering  model  we  are  using  to  compute  the  radiative  flux  in 
the  energy  equation  is  divided  into  two  separate  contributions.  First,  the 
direct  absorption  of  the  solar  energy  by  the  hydrogen/seedant  mixture  is 
computed  directly  by  solving  the  radiative  transfer  equation  along  discrete 
rays  of  concentrated  solar  beam.  The  only  approximation  here  is  a 
subdivision  of  the  broad-band  radiation  into  distinct  bands. 

The  second  contribution  to  the  radiation  flux  term  is  the  re-radiation 
by  the  hot  gases  and  the  confining  walls.  The  radiation  from  the  hot  gases 
is  being  modeled  by  a  differential  approximation  to  the  complete  equation 
of  transfer,  the  so-called  PI  approximation.  The  PI  equation  at  a 
particular  frequency  is  given  by, 


d_  3Gu  a_  x_  ^ 
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where  Gv  is  the  incident  radiation.  The  source  term  in  the  energy 
equation  is  then  computed  from  Gv  by, 

7'^  ■  kv  <4,Ibv  -  Gv)  I7) 

The  representation  of  the  frequency  dependence  requires  a  further 
approximation.  Standard  methods  for  incorporating  wavelength  variations 
include  the  Planck  and  Rosseland  mean  free  paths.  These  averages  have  the 
advantage  of  retaining  the  simplicity  of  the  grey  formulations,  but  neglect 
important  contributions  from  optical  windows.  The  alternative  is  to  employ 
a  simple  box  model  for  the  various  bands  in  the  system.  Here,  we  are  using 
a  band  model  along  with  the  PI  differential  approximation. 

The  spectral  information  from  the  detailed  description  of  the 
absorptivity  given  above  is  used  to  identify  the  appropriate  bands  and  an 
average  absorption  coefficient  is  used  in  each  band.  Preliminary  solar 
energy  absorption  calculations  in  a  flowing  gas  based  upon  a  four-band 
model  are  given  in  Figs.  17  to  19.  The  use  of  a  four-band  model  implies 
four  PI  equations  must  also  be  coupled  with  the  solar  energy  absorption  and 
the  energy  equation.  (For  these  initial  results,  the  flowfield  was  assumed 
to  be  one-dimensional  and  uniform,  pv-0,  pu-constant,  and  the  absorptivity 
variation  with  temperature  was  ignored.)  Figure  17  shows  heat  addition 
contours  for  an  energy  size  of  125N.  Figures  18  and  19  show  temperature 
contours  for  two  wall  temperature  cases.  Coupling  with  the  low  Mach  number 
equations  described  above  is  currently  being  initiated  along  with 
temperat ure-dependent  absorpti vi ty . 
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Fig.  1.  Convergence  of  transonic  time  dependent  algorithms  at 
various  Mach  number  conditions  for  a  two-dimensional 
flowfield.  \ 
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NUMBER  OF  ITERATION 


Fig  2.  Convergence  of  Euler  implicit  scheme  for  the  two  dimensional 
compressible  Euler  equations  using  preconditioning 

for  various  Mach  numbers,  CFL-5,  u>9»0.25.  Convergence 
plot  shown  is  for  the  continuity  and  x-momentum  equation. 
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Unmodified  Transonic 
Algorithm 


NUMBER  OF  ITERATION 


Convergence  rate  of  low  Mach  number  time  dependent 
algorithm  for  Mach  number  conditions  down  to  10  . 

For  comparison,  the  convergence  rate  of  the  original 
algorithm  is  shown  at  a  Mach  number  of  u.io. 


NUMBER  OF  ITERATION 


Convergence  rates  of  low  Mach  number  viscous  formulation 
for  various  Reynolds  numbers.  Present  results  are  for  a 
flow  Mach  number  of  10"4. 
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(b)  Velocity  Contours 

Re* 5000,  0*1  Re* 500,  0*1  Re*50,  0*4  Re=0.05,  0=4x10 

Flowfield  solutions  for  specified  heat  addition  in  a  duct  at  various  Reynolds  numbers. 

All  calculations  at  a  Mach  number  of  0.5  x  10~3. 


6.67 


Flowfield  solutions  for  flows  with  strong  wall  heat  transfer 
over  a  wide  range  of  Reynolds  numbers.  Flow  Mach  number,  0.5  x  10 


Mach  number -Reynolds  number  convergence  regime  for  regional 
time  dependent  algorithm  and  for  modified  low  Mach  number 
viscous  formulation  showing  approximate  regime  of  energy 
addition  problems  in  propulsion. 


Fig.  9.  Absorption  chamber  geometry  showing  stretched 
fluid  dynamics  grid  (51  x  31)  and  laser  ray 
comprising  a  doughnut-shaped  beam. 
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Temperature  Field  with  Velocity  Variation 
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Fig.  12.  Temperature  contours  in  laser  absorption  regime. 
Fixed  laser  power  and  varying  inlet  velocity. 


Temperature  Field  with  Chamber  Geometry 
Variation 
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Fig.  14.  Dimer  absorption  characteristics  of  sodium  at  500, 
2500,  and  5000K  based  on  four-band  model. 
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Fig.  15.  Dimer  absorption  characteristics  of  potassium 
at  2000  and  4000K  based  on  band  model . 


Fig.  16.  Photoionization  absorption  coefficient  of  sodi 
at  5000K. 
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Fig.  18.  Temperature  contours  for  non-gray  sodium  with 

absorptivity  independent  of  temperature.  T  =300K. 
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TEMPERATURE  CONTOURS  NON— GRAY  (CONST  K)  NA 


Fig.  19.  Temperature  contours  for  non-gray  sodium  gas  with 
absorptivity  independent  of  temperature;  T  =1000K. 
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System",  AIAA  Electric  Propulsion  Meeting,  New  Orleans,  LA,  November 
19,  1982. 

"Prediction  of  the  Flowfield  in  Laser  Propulsion  Devices",  AIAA  T8th 
Thermophysics  Conference,  Montreal,  Canada,  May  1983. 

"The  Effect  of  Strong  Heat  Addition  on  the  Convergence  of  Implicit 
Schemes",  AIAA  6th  CFD  Conference,  Danvers,  MA,  June  1983. 

"A  Two-Dimensional  Analysis  of  Laser  Heat  Addition  in  Converging 
Nozzles",  AIAA  Aerospace  Sciences  Meeting,  Reno,  NV,  January  1984. 

"High  Power  Nd-Glass  Laser  Instrument  for  Advanced  Propulsion  and 
Diagnostics",  AFOSR/AFRPL  Rocket  Propulsion  Research  Meeting, 
Lancaster,  CA,  March  12-15,  1984. 

"Stability  of  Absorption  Phenomena  in  Laser-Thermal  Propulsion 
Gasdynamic  Interaction",  AIAA  17th  Fluid  Dynamics,  PI asmadynami cs  S 
Lasers  Conference,  Snowmass,  CO,  June  25-27,  1984. 

"Analytical  Modeling  of  Strong  Radiation  Gasdynamic  Interaction" , 
AFOSR/AFRPL  Chemical  Rocket  Research  Meeting,  Lancaster,  CA,  March  20, 
1985. 

“Numerical  Solution  of  Strong  Radiation  Gasdynamic  Interactions",  AIAA 
18th  Fluid  Dynamics,  PI asmadynami  cs  &  Lasers  Meeting,  Cincinnati,  OH, 
July  16-18,  1985. 

"An  Implicit  Time-Dependent  Scheme  for  Low  Speed  Flow  with  Buoyancy", 

2 2nd  Annual  Meeting  Society  of  Engineering  Science,  Inc.,  University 
Park,  PA,  Oct.  7-9,  1985. 

"Computation  of  Compressible  Flows  at  Very  Low  Mach  Numbers",  AIAA 
Aerospace  Sciences  Meeting,  Reno,  NV,  Jan.  6-9,  1986. 
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"Computation  of  Low  Mach  Number  Flows  With  Buoyancy",  presented  at  the 
10th  International  Conference  on  Numerical  Methods  in  Fluid  Dynamics, 
Beijing,  China,  June  23-27,  1986. 

"Numerical  Solution  of  Strong  Radiation  Gas  Dynamic  Interactions", 
AFOSR/AFRPL  Chemical  Rocket  Research  Meeting,  Lancaster,  CA,  September 

1986. 

"Computation  of  Low  Speed  Viscous  Flows  with  Heat  Addition",  presented 
at  the  Sixth  IMACS  International  Symposium  on  Computer  Methods  for 
Partial  Differential  Equations,  Lehigh  University,  June  23-26,  1987, 
Bethlehem,  PA. 

"Coupling  Between  Concentrated  Solar  Radiation  and  Gasdynami cs" , 
AFOSR/AFRPL  Chemical  Rocket  Research  Meeting,  State  College,  PA,  June, 

1987. 
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V.  INTERACTIONS/ADVISORY  FUNCTIONS 


Member  of  Workshop  Panel,  "Concepts  and  Experiments",  NASA/Michigan 
State  Symposium  on  Advanced  Propulsion  Concepts  Using  Time-Varying 
Electromagnetic  Fields,  February  1982. 

"Aerospace  Propulsion  at  Penn  State",  presentation  to  General  Robert 
T.  Marsh,  USAF,  The  Pennsylvania  State  University,  University  Park, 
PA,  May  5,  1982. 

"Modeling  of  Flowfields  that  Interact  with  Radiation  Fields", 
AFRPL/UDRI  Solar  Plasma  Propulsion  Workshop,  Bergamo  Center,  Dayton, 
OH,  Jan.  21-22,  1986. 


VI.  NEW  DISCOVERIES 


